arXiv:gr-qc/0503108vl 25 Mar 2005 


Canonical quantization and black hole 
perturbations* 


Zoltan Perjes and Arpad Lukacs 


KFKI Research Institute for Particle and Nuclear Physics, 
Budapest 114, P.O. Box 49, H-1525 Hungary 
arpi@rrnki. kfki. hu 


February 7, 2008 


Abstract 

We examine the possibility of a constraint-free quantization of lin¬ 
earized gravity, based on the Teukolsky equation for black hole perturba¬ 
tions. 

We exhibit a simple quadratic (but complex) Lagrangian for the Teukol¬ 
sky equation, leading to the interpretation that the elementary excitations 
(gravitons bound to the Kerr black hole) are unstable. 
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1 Introduction 

Quite a number of attempts have been made and are being made to quantize 
Einstein’s gravitation theory. One of the earliest attempts to use canonical 
quantization of linearized gravity goes back to Gupta [Q. The canonical ap¬ 
proach to quantum gravity using the 3+1 decomposition has been worked out 
by Arnowitt, Deser and Misner 2]. A different idea has been put forward by 
Penrose [5] introducing twistors. Ashtekar 0j introduced spinor variables to 
handle the constraints. Nowadays loop quantum gravity is being actively pur¬ 
sued, see the review by Rovelli 0. In spite of the remarkable progress achieved 
using these beautiful constructions, there remains a lot to be done. 

Undoubtedly most of the difficulties in all canonical approaches to quantum 
gravity come from the constraints. Another problem, closely related to the 
treatment of constraints is the choice of the canonical coordinates. A good choice 
of coordinates leads to less constraint equations, and thus it makes quantization 
easier. 

In general relativity, the physical state of space-time is described by a sym¬ 
metric, second rank tensor that can be parametrized by ten components. Due 
to the constraints, we should introduce less than ten independent canonical 
coordinates, and their conjugate momenta. 

*Talk presented by Zoltan Perjes 


1 




The well known Teukolsky equation (Bjj, which is a linear PDE, describes the 
perturbations of a Kerr space-time with a single scalar quantity, ip, one of the 
Newnran-Penrose spin coefficients. 

Chrzanowski, Misner fj and Ori {5, have shown the way to reconstruct the 
metric perturbations from the solution of the Teukolsky equation by acting on 
it with a differential operator. 

This makes the quantity ip a promising candidate as a canonical coordinate 
in a linearized quantum gravity, provided a Hamiltonian formulation of the 
Teukolsky equation IB can be given. This problem will be examined in the 
next section. 

One would think, that the easiest way to obtain a canonical formalism based 
on the Teukolsky equation is to reconstruct the perturbations of the metric ten¬ 
sor, and subsequently the quadratic part of the corresponding perturbed La- 
grangian. This could then be taken as a starting point for a canonical quantiza¬ 
tion. However, the formulae for the metric perturbations are very complicated, 
and contain fourth order derivatives of ip. 

One could just as well try to find a Lagrangian directly for Teukolsky’s 
equation in itself, and use that for quantization. Unfortunately, as we will 
see later, there is no quadratic, first order real function, which gives IB as an 
Euler-Lagrange equation. 

We have identified a complex Lagrangian for the Teukolsky equation. A 
possible solution to quantize the Teukolsky equation could be considering the 
complex part of the Lagrangian as an interaction term. The interaction-free part 
can be then simply quantized, and it’s quanta may be identified as gravitons. 
As a result of the perturbations (the complex part of the Lagrangian), these 
excitations decay. 


2 Is a Lagrangian formulation of Teukolsky’s equa¬ 
tion possible? 


Some basic features of the Teukolsky equation can be found in the Appendix. 
Here we recall that this equation governs the dynamics of one of the Riemann 
tensor’s tetrad components in a linearized approximation over a Kerr back¬ 
ground. 

Using the notations of the Appendix, Teukolsky’s equation (Hj can be written 
as 
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(1) 

where T is a source term, calculated from the energy-momentum tensor of other 
fields. In the following, we concentrate on the free gravity case (T = 0). 

Over a Kerr geometry, this equation describes massless scalar fields (s = 0), 
electromagnetic perturbations (s = ±1) and gravity (s = ±2). 
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In this section, we examine, whether it is possible to obtain equation 0 as 
an Euler Lagrange equation from a real Lagrangian of the form 

C = ag lk di(bp)d k (bip) + cpip, (2) 


which is the most general first order quadratic local Lagrangian of a scalar field. 
The motivation to try to find such a Lagrangian (eg. local, not containing higher 
derivatives) is the possibility to quantize it easily. 

In order for C be real, a and c have to be real, and b complex arbitrary 
functions of the coordinates. From the above function, the Euler-Lagrangian 
equation is given by the formula 


Q dy/^qC _ d^f^gC 
ddiip dp 


(3) 


where 

\f—g = (r 2 + a 2 cos 2 i9) sind, 

is the determinant of the metric tensor. 

Calculating the above equation, and taking the coefficients of p’s derivatives 
in 0 and 0 yields equations for a, b, c and their derivatives. When trying 
to solve these equations, one finds, that the integrability conditions (eg. the 
r derivative of da/dt should be equal to the t derivative of da/dr) cannot be 
fulfilled. This clearly shows that no real Lagrangian of the form 0 can lead to 
Teukolsky’s equation. 

While this work was in progress, we have learnt from the summary of Stephen 
Anco’s talk at the Montreal Workshop on the Interaction of Gravity with Exter¬ 
nal Fields 0, that Teukolsky’s equation has no Lagrange function. According 
to Anco 0, this is because the non self-adjointness of 0 , which means that 


J pDp^ J pDp, 

where D denotes the differential operator in 0 . 

To show that the equation is indeed not self-adjoint, one should take the 
terms in f pDip and do the partial integrations to calculate its adjoint, and 
compare it with the complex conjugate of f ipDp. 

In the first term, self-adjointness is obvious: the coefficient of d 2 /dt 2 does 
not depend on t. Similar is the case with the second and the third term. 
However, the non-self-adjointness can be seen on the fourth term: 
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where substituting A, and expanding the derivatives yields in the s = 2 case 
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when doing the partial integrations. One can see, that exchanging the roles of 
i[) and <f> in the original integral, and conjugating does not give such terms, and 
no other term can cancel the r-derivatives in the above. 

Bini, Cherubini and Jantzen m gave another form of the above equation, 
closely resembling an ordinary Klein-Gordon equation. They have proven, that 
the Teukolsky equation is a generalization of the Klein-Gordon equation for 
arbitrary spin weight s. Examining the following form of the equation, it is 
easy to point out the cause of the non self-adjointness: 

[(V 4 + sT i )(Vi + sTi) - As 2 ^] i/> = 4ttT, (6) 


where T is a four-vector with components 
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It is easy to show that 
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Now we can see, that the equation is non self-adjoint, because the p- 
components are not purely imaginary, unlike the Klein-Gordon case, where 
Tj = iAi, where A, is the electromagnetic four-potential. These complex terms 
make the equation absorptive. 


3 Perturbative quantization of the Teukolsky equa¬ 
tion 

In the previous section we have seen that a direct canonical quantization of 
the Teukolsky equation is not possible because of the absorptive terms in the 
Lagrangian. On the other hand, the Bini-form © of the Teukolsky equation 
© is in close resemblance to the Klein-Gordon equation. This suggests the 
following solution: we separate the absorptive part of the equation to treat it as 
a perturbation, and this way we can construct a Lagrangian for the self-adjoint 
part. This will be real, and so the Fock space construction is well defined, a 
particle interpretation can be given. 

Introducing the notation of 


Ai = 9f(T0 

Bt = R(T0, 
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the Bini form of the equation can be written as 

[(V* + s(B 1 + *A i ))(V i + a(Bi + iAJ) - 4 s 2 ^] = 4 ttT, (9) 

and if we keep only the -AjS, we get an equation which is formally a Klein-Gordon 
equation in an external electromagnetic field. The (complex) Lagrangian for this 
equation is: 

C = {V i -s(B i + iA i )}^{\/ i + s(B i + iA i )}^ + 4sH^, (10) 

which can be written as a sum of an unperturbed part, and a perturbation term: 

C = Ckg + Ci, 


where 

Ckg = (V* + isA l )ij)(y i + isAi)%'> 
is the Klein-Gordon Lagrangian, and 

C\ = — B l il)Viip) + 2isB’ l Ai'ijj'il) + 4 zs 2 4'^i/’V’ 

is the perturbation term. 

The canonical momenta are given then by the well-known formula 

n = ~^= = (<9° + isA°)ip 
dd 0 $ 

n = WT = ( 5 ° - iaA °ffl’ 

od 0 ip 

and the Hamiltonian can be expressed as 

Ho = n<9oV> + n<9 0 ? - £ 0 = nn + (v Q + iA a ) 4 >{v a + iA a )i/j. 


( 11 ) 

( 12 ) 


(13) 


(14) 


In the following discussion Ho will be taken as the ’’unperturbed” Hamilto¬ 
nian density for gravitons. In the next section, we examine the ’’unperturbed” 
equation obtained by taking the variation of the Lagrangian Ckg- 


4 The ” unperturbed” equation. Separation of 
variables 

The Euler-Lagrange equations can be obtained by using either 0 or the well 
known tensorial - and thus simpler - formula 

_dCo =0 

dViij) di/j 

The resulting equation, as mentioned above, is formally a Klein-Gordon equa¬ 
tion with an external electromagnetic four-potential: 

(Vi -isAi)(V l -isAJip = 0 (15) 
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In the following paragraphs, we examine this equation in detail. An explicit 
computation yields 


((( 2M - r ) r - « 2 ) + 0 - 4 J0 aMr ~ { ° 2 ~ 2Mr+r2) 0 
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v ' at 2 \ at ox) ap> J ' J 

( 16 ) 

It can be shown that the variables in the above equation are separable. Let us 
look for solutions in the form 

r, d, ip) = e~ iut R(r)T(d)e irnv . (17) 

The Euler-Lagrange equations can be normalized so, that the coefficient of 

o2 I __ 

will be 1, similarly to the original Teukolsky equation ifljl. 

After substituting the ansatz o into the above equation, it can be seen, 
that none of the terms contains mixed derivatives, and the angle dependence 
can be factorized in the coefficients of the radial derivatives, and similarly the 
r-dependence can be factorized in the coefficients of the ^-derivatives. 

The above properties of the equation make the separation of the variables 
possible. Straightforward but tedious calculation yields for the radial equation 

(a 4 - 4 a 2 Mr + 2a 2 r 2 + 4AfV - 4Mr 3 + r 4 )R"(r) 

(IQ) 

+2(— a 2 M + a 2 r + 2 M 2 r — 3Afr 2 + r 3 ))R’(r) = A r R(r), 


where 

X r = a 4 w 2 — Aa 2 Mw 2 r — a 2 m 2 — a 2 s 2 — AaMmwr + 2 Mrs 2 — w 2 r 4 — r 2 s 2 — k , 

in which k is the separation constant. 

Similarly the angular equation is 

(cos 2 d — 1 )T''(d) — cost?sin dT'(d) = A $T(d ,) (19) 


where 

A,j = ((cos !?aw — s ) cos 2 daw + (aw + m)s ) cos d + m 2 + s 2 — arw 2 + k. 

Introducing x = cos d as a new variable in the angular equation we get 

(— x 4 + 2a; 2 - 1 )T"(x) + 2x(l - x 2 )T'(x) = A $T(t). (20) 

Both equations are boundary value problems for second order linear differ¬ 
ential operators, with regular boundary conditions at the horizon and at r = oo. 
Unfortunately we have not succeeded to express the solutions of equations m 
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and CD using elementary functions. Let ip k denote the kth eigenfunction of 
m, then the general solution can be written as 


^=y^Qfc^fc- (21) 

k 

Now the canonical quantization of the unperturbed problem is completely 
straightforward. The modes m of equation m can be interpreted as gravitons 
bound to the black hole. Defining the corresponding a k and aj^ absorption and 
emission operators for each mode, the Hamiltonian can be expressed as 

n = J2 E k a\a k . (22) 

k 

The E k energy levels can be determined numerically, as well as the eigenmodes. 

Taking into account the complex part of the Hamiltonian (1121) of the full 
Teukolsky equation, the ’’gravitons” of (EU become unstable. This may not 
come as a complete surprise, as we know that the Kerr black hole admits only 
quasi normal modes. 

The decay constants for these excitations can be calculated using the ma¬ 
trix elements of the perturbation part of the Lagrangian, similarly to the way 
decaying states are described with an imaginary energy contribution in nuclear 
physics. 


Appendix 


In this paper, we quantize the perturbations of a Kerr black hole. We use the 
Boyer-Lindquist coordinate system: the coordinates are the time, f, the radius, 
r, and the usual spherical angles, d and ip. Using these coordinates the line 
element can be expressed as 
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where we used 

and 
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for simplicity. Here M denotes the mass of the black hole, and 

J 
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is the fraction of its angular momentum (J) and mass. It’s an important prop¬ 
erty of these coordinates when deriving the Teukolsky equation, that these han¬ 
dle the two double principal directions of the space-time symmetrically. 
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To derive the Teukolsky equation, the Newman- Penrose null tetrad method 
is utilized. When using Boyer-Lindquist coordinates, a usual choice for the null 
tetrad is the Kinnersley tetrad: 
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Using these, the non-vanishing spin coefficients are 
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Here we recall the basic properties of the Teukolsky equation 0. This equa¬ 
tion is obtained using the Newman-Penrose tetrad method. This calculation 
can in fact be done for any space-time in the D Petrov class (ie. two pairs of 
degenerate principal directions). 

The perturbations of a Kerr space-time are described by the Teukolsky equa¬ 
tion 
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then. The calculation of T from the energy-momentum tensor of other fields is 
not included here, but it can be found in Teukolsky’s original paper [5j. 

In this equation, the variables can be separated, using an ansatz 


i/> = e~ iut e im,fi S(e)R{r) 


for the dependent variable. The radial equation is 
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and the angular equation is 
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The boundary conditions demand %l> being regular at d = 0,7r. The solutions 
are examined in a paper by Mano, Suzuki and Takasugi Ei. 

According to Bini, Cherubini, Jantzen and Ruffini cm, the above equation 
is the direct generalization of the Klein-Gordon equation for arbitrary s spin 
weight. In their paper, they have brought the equation into a form closely 
resembling the Klein-Gordon equation 


[(V 1 + sT l )(V i + sTi ) - 4 S 2 4^] if, = 4t tT, 
where T is the following four-vector: 
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